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1 Introduction 

Interface states have been extensively studied in rigorous statistical mechanics since Dobrushin's 
seminal paper on interfaces Gibbs states in the three-dimensional Ising model W. Interfaces and 
surfaces are everywhere in the physical world, and to know their properties is important in un- 
derstanding phenomena as diverse as droplet evaporation, crystal growth, catalysis, dynamics of 
magnets, etc.. Until recently almost all rigorous work was restricted to classical models, while some 
of the phenomena are certainly more adequately described by quantum mechanical models. In the 
past few years, however, several interesting results have been obtained for the quantum Heisenberg 
and the Falicov-Kimball model. The Falicov-Kimball model is discussed in several other contribu- 
tions in the this volume [§,[§-) Therefore, we restrict our attention here to the Heisenberg model 
and quantum spin systems in general. We will also introduce a quantum solid-on-solid (QSOS) model 
for the 111 interface of the Heisenberg XXZ ferromagnet. 



2 Interface ground states of the XXZ ferromagnets 

A quantum spin model is defined by specifying a dynamics on an algebra of quasi-local observables. 
The local structure is given by the finite subsets of the d— dimensional latttice 7L d . With each site 
x G Z d , there is associated a copy A x of the n x n matrices with complex entries M n (C). A x is the 
algebra of observables at the site x. For every finite subset A C Z d , the observables in the volume 
A are given by 

This is a finite-dimensional C*-algebra, and if Ao C A, we have the natural embedding 

'Copyright © 1998 by the authors. Reproduction of this article by any means is permitted for non-commercial 
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The algebra of local observables is then defined by 

■Aloe = (J A A 

AcZ d 

where the union is over all finite subsets of Z d . Its completion is the C* — algebra of quasi-local 
observables: 

A — Aloe 

We will also need the translation automorphsims on A, denoted by t x ,x G Z d , canonically 
mapping Aa into ^a+i- 

The dynamics is determined by a family of local Hamiltonians. For simplicity we will only discuss 
models with translation invariant finite range interactions. I.e., let h — h* G Aa 01 for some finite 
set Ao, and define the local Hamiltonians by 

H A = T *W 

i:Ao+iCA 

The generator of the dynamics is the unique closed extension of the derivation 

8(A) := lim \H A , A], for all A G A\ oc 

A\l d 

This generator can be exponentiated to obtain a strongly continuous one-parameter group of C*- 
automorphisms {a t }t£R, 

a t (A) := e it5 {A), A e A 

The standard proof can be found in the books by Bratteli and Robinson j|] or Simon Q . 
A state of the quantum spin system is a linear functional uj on A with the properties: 

uj(A*A) > 0, for all A G A, lu(1) = 1 

The ground states of a model, by definition, are the solutions of the following set of inequalities: 

lu{A*8{A)) = lim u(A*[H A ,A]) > 0, for all A G Ai oc 

Atoo 

These inequalities express that all local excitations, i.e., created by a local observable A, raise the 
energy. This is also called Local Stability and we will often refer to this set of inequalities as LS. 
For translation invariant states, i.e., states satisfying 

oj o t x — uj, for all x G Z d 

the ground states are exactly the states that minimize the energy per site: 

uj(h) = mf{ri(h) \ r\ translation invariant state on .4} 

In general, LS has non-translation invariant solutions, e.g., describing domain walls or interfaces. 
In this paper, our main interest are these non-translation invariant solutions. In one dimension 
they are often called kink (and antikink) states, or soliton states. Solutions of LS are commonly 
constructed as limits as A — > Z d , of ground states of the finite- volume Hamiltonians plus boundary 
terms. 

One-dimensional quantum spin models, also called spin chains, are of particular interest for 
several reasons. First of all it is in one dimension that the most detailed rigorous analysis is 
possible. This was the main motivation for 0. More recently it has also become possible to realize 
quantum spin chains experimentally and to compare theory and experiment in surprising detail. It 



Interface states of quantum spin systems 



3 



is also interesting that there is a special but not so small class of one-dimensional models for which 
the exact ground states can be given explicitly [|[ ||, |l0[ [n]]. And finally there are the integrable 
quantum spin chains and the rich mathematical structures they exhibit jl2j . 

None of the works cited above, however, deal with non-periodic ground states. Consequently, 
the problem of characterizing the complete set of solutions of LS, for any model, was, until recently, 
never dealt with. 

The most detailed results have been obtained for the ferromagnetic XXZ Heisenberg model. For 
x G Z d , let S x , i = 1,2,3, be the standard spin-5* matrices, generating a n = 2S + 1-dimensional 
irreducible unitary representation of SU(2), with S — 1/2, 1, 3/2, . . .. The local Hamiltonians of the 
ferromagnetic XXZ Heisenberg model are given by 

Ha = - ^(SlS'y + S 2 x S 2 y ) + S 3 x S 3 y (2.1) 

\x — y| — 1 

with A > 1. A = +oo is the Ising model. A = 1 is the isotropic model, also called the XXX model. 

First, we discuss the complete set of ground states, in the sense of LS for the XXZ ferromagnetic 
chain, i.e., d = 1. This class of models depends on two parameters: the dimension of the spin 
matrices (25* + 1, S = 1/2, 1, 3/2, . . .), and the anisotropy A > 1. 

Consider finite volumes of the form A = [a, b] C Z, and denote by <9A = [a, a + r] U [b — r, b], the 
boundary of A, where r > is a suitably chosen integer. 

Solutions of the ground state inequalities can be constructed by adding suitable boundary terms 
&a € Aqa to the finite-volume Hamiltonians, and taking limits A | Z of finite volume ground states 
of the form 

ca(A) = ^M {0Ta]1 AeA A 

{VA,Va) 

where ipA £ 0igA^ 2S+1 i ^ s an eigenvector belonging to the smallest eigenvalue of H\ + b\. 
For the XXZ chains it suffices to take r = 1 and boundary terms of the form 

b [a , b] = B(S 3 b+1 - S 3 ^) (2.2) 

For convenience we will consider the chain on intervals of the form [-L, L], including the boundary 
points (—a = b = L — 1). 

Two translation invariant solutions are trivial to find and have been well-known for a long time: 
the unique states a>j and u>± determined by 

w T (££)= S, forallxeZ 
Ul(Sl) = -S, forallxeZ 



By taking B := B(A) ± Sy/l — 1/A 2 , Alcaraz, Salinas, and Wreszinski |u|, and independently, 
Gottstein and Werner fl4| , found non-translation invariant solutions. A first parameterization of 
these kink states is by the third component of the total spin, i.e., the eigenvalue M of J2x=-l ^ii 
which commutes with the Hamiltonian. The possible values of M are —S(2L + 1), —S(2L + 1) + 
1, . . . , S(2L + 1). For each value of M there is exactly one ground state -0m- It is useful to define a 
generating function, tp(z), zeC, 

ip(z) = ^ a m z M ipM (2.3) 

M 

where the <zm are suitable constants. We refer to |l^, O for more details. The important fact is 
that the au can be chosen such that (p(z) becomes a product: 

L 

¥>(*) = 6$ X(zq~ x ) (2-4) 
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with 



2S+1 



s 



m=-S 



z q 



x(S- 



(2S) 



y (S - m)\{S + to)! 1 ' 

where S 3 \m) = m \m), and q G (0, 1), such that A = (q + q~ 1 )/2. Note that, in the limit L —> oo, 
a translation by one lattice unit transforms tp(z) into ip(zq). Therefore, we can use tp(z) as a 
parameterization of the kink states by translations. It is then a trivial computation to get the 
magnetization profile in the states ip(z). E.g., for S — 1/2, one finds 



Wkink(S x ) 
Wantikink(S'l) 



-tanh((x-a)/£) 
--tanh((a:-a)/f) 



where £ and a are simple functions of q and z. 

The concept of zero- energy states plays an important role in the characterization of the complete 
set of solution of LS. So, we introduce it here, before we continue our discussion of the XXZ model. 

Suppose there exists < h G ^Iaq , such that 



Then, if a state ui satisfies 



H A + b A := H A = T *Ch) 

x:A +xCA 



w(r x (/i)) = 0, for all x G Z, 



u> is called a zero-energy state. It is easy to show that any zero-energy state satisfied LS. 

Gottstein and Werner obtained the complete set of zero-energy states of the anisotropic XXZ 
chain A > 1 with the particular boundary terms given in ( |2.2j ). They proved that any pure zero- 
energy state is either one of the two translation invariant ones, or it is a member of a set of mutually 
equivalent kink states or a set of mutually equivalent antikink states. A state u> is called pure, iff 
for any two states wi, u>2, and t £ (0, 1), one has 

UJ = tu>\ + (1 — t)u>2 LOi = U>2 

It is obvious that the solution set of LS is convex, and one can prove it is a face. The same holds 
for the set of zero-energy states of a fixed h. Therefore, finding the pure solutions is enough. 

To which class u> belongs is determined by the limits a, G {±5}, i.e., its asymptotic behavior: 

a := lim ui(SZ), (3 := lim u>(Sz.) 

x — > — oo x — > + oo 

The following table summarizes the four types of zero-energy ground states as parametrized by a 
and P (A > 1): 



a (3 


type 


dominating configuration 


++ 


up 






down 




-+ 


kink 


1 • • • 1 I I I I ••• T T T T T T •••) 


+- 


antikink 


1 • • • T T T T T ••• I I I I I I •••> 



We refer to @ for more details on the kink and antikink states. The case A = 1 will be 
discussed further on. 

Let us now return to the question of characterizing all solutions of Local Stability. The first 



result is due to Matsui 15 
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Theorem 2.1 (Matsui, 1996) For the ferromagnetic XXZ chain, with A > 1, and S = 1/2, the 
translation invariant states cj-j- and u>i, together with the kink and antikink states described in JT^ / 
and discussed above, are the full set of pure solutions of LS. 

Note that this theorem does not say anything about the isotropic model (A = 1). The isotropic 
ferromagnetic chains have an infinite family of translation invariant ground states, given by the state 
lui and all rotations of it. As expected, this breaking of the rotation symmetry is accompanied by 
a gapless excitation spectrum. The proof of the above theorem relies on the existence of gap, and, 
therefore, does not extend to the isotropic case. Recently, Koma and Nachtergaele obtained a proof 
of the complete set of ground state for all values of S and A > 1 Jl(| . 

Theorem 2.2 (Koma and Nachtergaele, 1997) For the ferromagnetic XXZ chain, with A > 1, 
and S > 1/2, the translation invariant states uj^ and £Jj_, together with the kink and antikink states 
generalized to arbitrary S, as in JZj}/ exhaust the set of pure solutions of LS. If A = 1, all solutions 
are translation invariant. 

So, in the isotropic case there are no kink-type ground states. One expects that the same is true 
for models with a unique translation invariant (or periodic) ground state. This has not yet been 
proved in general, but Matsui [|l7) obtained a quite general result in the case of a unique zero-energy 
state, which we explain next. 

Consider an arbitrary spin chain with a translation invariant nearest neighbor interaction, i.e., 
h e A[ ,i] : 

6-1 
x—a 

Assume that h > and that there is a unique translation invariant state u> such that 

u{h x ,x+i) = 0, for all igZ 

One would expect that in this case uj is the unique solution of LS. Matsui's result, the only result so 
far, requires an additional assumption on the set of zero-energy states of the the half-infinite chains, 
i.e., states r/ on «4[i.+oo) (•A(- O o,o]i respectively) such that 

ili^x.x+i) = 0, for all x > (x < 0, resp.). 



Theorem 2.3 (Matsui, 1997) For a spin chain as described above, if to is the unique zero-energy 
ground state and if all zero-energy states of the left and right half-infinite chains are quasi- equivalent, 
then lo is the unique solution of LS. 

In particular this implies that the AKLT spin 1 chain introduced in Q has a unique ground state 
in the sense of LS. 

In higher dimensions there also exist non-translation invariant ground states. In finite volume, 
this was already noted in fl3fl . Koma and Nachtergaele proved that ground states with a rigid 
interface in the diagonal (i.e., 11 • ■ • 1) direction exist in all dimensions. They also proved that the 
excitation spectrum above the ground state with a 11 interface in two dimensions is gapless [ff8[ . 
Matsui generalized this to arbitrary dimensions d > 2 jl9| . We now give an easy construction of a 
111 interface state for the XXZ ferromagnet. 



It follows from ( |2.4| ) that, for all z S C, 

h x (z) ® x(z<i~ l ) = o, 



where 

h := -(S 1 <g> S 1 + S 2 ® S 2 ) - A(S 2 <g> S 3 - S 2 ) + B(A)(S 3 <g> 1 - 1 <g> S* 3 ) (2.5) 
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Figure 2.1: The decomposition of the square lattice into coupled zig-zag chains. The 
full lines are the bonds in the chains. The dashed lines are bonds carrying the interchain 
couplings. The triple dotted line indicates the position of the diagonal interface. All vertical 
bonds are oriented upwards, and all horizontal bonds are oriented to the right. 



is the positive definite interaction for which the kink states are the zero-energy ground states. Let 
h x —>y denote a copy of h acting on the Hilbert space Ti.{ x , y } of the sites x and y in this order. It 
follows immediately that, for any A C Z d , x, y £ A, such that \x — y\ — 1, we have 

h x ^ y Vt{z) = 

with 

:=(g)x(^<r H ) 

tiGA 

Hence, 

H A n(z) = h x ^ y fl(z) = (2.6) 

x lH eA,|an|<| H 
|«-t/|=l 

With the observation that every interior point in A is the starting point (a;) and the end point (y) 
of an equal number of bonds, it is easy to see that H\ is the XXZ Hamiltonian with a boundary 
field. We conclude that fl(z) is a zero-energy ground state for the XXZ model on Z d for any d. 

For concreteness, let us continue the discussion in the two-dimensional case. Generalization to 
arbitrary dimensions is straightforward. The two-dimensional system can be decomposed into a 
collection of coupled one-dimensional XXZ models defined the zig-zag chains shown as full lines in 



Figure |2.1| . The bonds coupling two zig-zag chains together are indicated by dashed lines. As f2(z) is 
a zero-energy ground state of the two-dimensional system, its restriction to a zig-zag line must be a 
zero-energy ground state of the one-dimensional XXZ model. Note that the boundary terms coincide 
with the ones used before to obtain the kink ground states if we interpret the orientations of the 
bonds in two dimensions as left-to-right in one dimension. Further inspection shows that the kink 
states on the zig-zag lines are aligned such that they are all centered on a line in the main diagonal 
direction (11) of the two-dimensional square lattice. Thus, £l(z) describes a 11 ■ • -1-interface. 
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Theorem 2.4 (Koma and Nachtergaele, unpublished; Matsui 1997) If d > 2, S > 1/2, 
and A > 1, the family of product states O(z), z£C, normalized and extended to infinite volume, are 
all the pure zero-energy ground states of the interaction defined in The excitation spectrum 

above these ground states is gapless. 

We already noticed that \z\ determines the position of the kinks, i.e., a change in \z\ corresponds 
to a translation of the interface. It is not hard to see that the phase of z can be changed by applying 
another symmetry of the Hamiltonian: 

e ies ^Q(z) = n{e l6 z) 

Rotations about the z— axis are commute with the Hamiltonian, but the states ft (z) are not invariant 
under it. Therefore, we have identified a broken continuous symmetry that goes along with the 
breaking of translation invariance in interface states. This is the basis of Matsui's general proof of 
the existence of gapless excitations above these interface states of the XXZ model, along the lines 



of 1 20 



For sufficiently large A, and d > 3, interface ground states (as well as low-temperature states) 
with the interface along the main coordinate directions (100), have also been shown to exist |2l], p2[ . 
In this case, however, there is a gap above the ground state. Another difference is that interfaces 
in the 100 directions can be understood as perturbations of the same interface states for the Ising 
model, which describes the limit A — > oo. The Ising model does not have a rigid 111 interface. It 
is easy to see that the ground states with 111 boundary conditions is infinitely degenerate. See [^3| 
for a discussion. For finite A, however, the XY term in the XXZ Hamiltonian produces a ground 
state selection effect, i.e., lifts the degeneracy. 



3 A quantum SOS model 

The next natural question concerns the existence of a rigid 111 interface at non-zero temperature. 
The structure of a 111 interface in a lattice model on Z 3 is already quite complicated. In the case 
of the XXZ model, the presence of gapless excitations in the interface ground states, is additional 
difficulty. In the past, interesting progress with interface problems has been made by introducing 
models where the interface is described by a height function. This eliminates most of the complicated 
geometry of the interfaces, which usually does not substantially affect the most relevant physical 
properties. Such models are called Solid-On-Solid (SOS) models. In this section, we introduce 
an quantum SOS (QSOS) model for the 111 interface of the anisotropic Heisenberg ferromagnet. 
As this model is a quantum system, the most straightforward and reliable way of obtaining an 
approximation is by restricting the Hilbcrt space to a suitable subspace. We would also like to 
retain as much as possible particular features, as, e.g., the gapless excitations of the interface. 

The restricted Hilbert space for the QSOS model should at least contain the ground state in- 
terface, and capture the most important excitations. In analogy with the classical SOS models 
we would also like to have the property that each one-dimensional subsystem defined on a line in 
the 111 direction, i.e., perpendicular to the interface, is found in one of its ground states. In the 
classical case this would mean that any such line crosses the interface exactly once in any SOS 
configuration. This implies that the interface can be described by a height function. In our case the 
one-dimensional subsystems are naturally defined on the zig-zag lines drawn in full in Figure 2.1. 
This idea can easily be generalized to 3 and more dimensions. For each one-dimensional the natural 
Hilbert space is then given by 

Hhz = closed linear span of {<f(z) \ z G C} (3.1) 

By going back to the basis of kink states with fixed magnetization M, one can easily show that this 
Hilbert space is also given by 

H zz = closed linear span of {<p(q n ) | n £ Z} = l 2 {2) 
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The states </?(<?"), n £ Z, are an non-orthogonal basis for 7i zz , and the index n can be interpreted 
as the height function, i.e., the position of the interface. There is a unitary operator 5* on TL ZZ such 
that 

^(0 = ^ ( " +1) ) (3-2) 

The Hilbert space for the QSOS model is then a tensor product of copies of H zz , one for each zig-zag 
line: 

HqsoS = H zz . u 

uer 

where u G T label the zig-zag lines. In the case d — 2, T is a simple one-dimensional lattice. For 
d = 3, r is a two-dimensional triangular lattice. 

The next step in the definition of the QSOS model is the calculation of the Hamiltonian. Let 
Pqsos denote the orthogonal projection onto Hqsos considered as a subspace of the original Hilbert 
space of the d— dimensional XXZ model. Then, the Hamiltonian of the QSOS model is given by 

Hqsos = PqsosHxxzPqsos 

Some symmetry properties follow immediately from this definition. The first is invariance under a 
global shift of the interface. This is expressed by the fact that 

[Hqsos,® S] = 
uer 

where S is the unitary operator defined by ( |3.2| ) . Because of this invariance, one needs to consider 
suitable boundary conditions when taking the thermodynamic limit in the directions parallel to the 
interface. One can also implement this by defining the infinite-volume model on incomplete tensor 
product in the sense of Guichardet [Q. Furthermore, from the definition of H zz it is clear that 
Hqsos is invariant under global rotations about the z— axis. Therefore, the rotation invariance of 
Hxxz is inherited by Hqsos- Finally, the translation invariance of the original Hamiltonian, which 
is also at the origin of the invariance under S (translations perpendicular to the interface), leads to 
periodicity or translation invariance of the QSOS model over T. 

There are different possibilities for the concrete representation of Hqsos ■ in order to represent 
this Hamiltonian as a matrix in an orthonormal basis one needs to correct for the non-orthogonality 
of the natural basis {ip(q n ) \ n G Z} of the TC ZZ . Although there are several examples in the 
physics literature where in the definiton of projected models non-orthogonality of a natural set of 
basis vectors is ignored, we think it is important to take it into account here. We will elaborate 
on this in a future publication. We conjecture that the QSOS model, with d > 3, has a rigid 
phase at sufficiently low temperatures. We will report on partial results on this problem in a future 
publication ]25| . 
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